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The bosonic IIB matrix model, which contains the bosonic part of the IIB matrix model conjectured to be a 
non-perturbative definition of the type IIB superstring theory, is studied using a numerical method. The large N 
scaling behavior of the model is shown performing a Monte Carlo simulation. The expectation value of the Wilson 
loop operator is measured and the string tension is estimated. From the numerical results, the prescription of the 
double scaling limit is obtained. 



1. Introduction 

Recently, several models have been proposed as 
non-perturbative formulations of string theory jl]- 
[|. The IIB matrix model 0, which is zero- volume 
limit of the ten-dimensional large N supersym- 
metric Yang-Mills theory, has been proposed as a 
constructive definition of the type IIB superstring 
theory. It is defined by the following action, 

s = ~7 tr (i [Al ' Av]2 + l^ A ^) ' W 

where and i\) are N x N traceless Hermitian 
matrices. The interesting feature of this model is 
that the space-time coordinates are considered as 
the eigenvalues of these matrices. Then, we ex- 
pect that the fundamental issues including the di- 
mensionality and the quantum gravity can be un- 
derstood by studying the dynamics of the model. 

To take the continuum limit (N — * oo) for the 
IIB Matrix model, a sensible double scaling limit 
should be determined dynamically for the scaling 
property of the two important quantities of the 
model, the string scale (a 1 ) and the string cou- 
pling constant (sw)f|]. 

a' ~ g 2 N a+b = g 2 N^ ~ Constant, 
9str ~ N a - b , (2) 
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where one have a restriction, a = b for the finite 
value of the string coupling constant {g s tr)- 

For studying the dynamical aspects of the 
model, some Monte Carlo simulations have been 
performed. In Ref. || , the bosonic model has been 
studied with a analytical method as well as a nu- 
merical one. The numerical results support the 
l/D expansion as an effective tool detecting the 
large N scaling behavior of the model. The lead- 
ing term of the 1 /D expansion at D > 3 suggests 
that the exponent (7) takes a value of 1. 

In analogy with the two-dimensional model, the 
Eguchi-Kawai model, we study the scaling prop- 
erty of the Wilson loop in ten-dimension. The 
area law of the Wilson loop operator has been 
found in the four-dimensional model in Ref. [Q. 
Then we consider that the scaling property of 
the string tension is estimated in ten-dimensional 
model with investigating the scaling property of 
the Wilson loop in ten-dimension. 

In this article, we present the data which show 
the existence of the double scaling limit in the 
model and discuss our numerical results. 

2. Monte Carlo simulation and scaling be- 
havior 

For numerical simulation of the model, we con- 
sider a bosonic version of the IIB matrix model, 
namely a model which can be obtained from the 



IIB matrix model by dropping the fermions by 
hand[^| . In the IIB matrix model, supersymmerty 
is expected to be very important as is discussed in 
many places of the literature j|,0-f§- However, it 
is difficult to calculate with numerical method. In 
ref. II , numerical study of a supersymmetric ver- 
sion of the IIB matrix model is discussed. In this 
article, we study the most simplest model. The 
bosonic model of the IIB matrix model is given 
by 

Sbosonic — otr^uj A v ] [A^ , A v \ , (3) 
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where are N x N Hermitian matrices rep- 
resenting the ten-dimensional gauge fields. The 
coupling constant (g) is nothing but a scale pa- 
rameter and is absorbed with the rescaling of the 
gauge field (A M ) as A^ -> ^=^V 

we consider the partition function of bosonic 
model in ten-dimension, 



Z = dAe 



-4 r tr[A )J ,A„][A' i ,A'' 



(4) 



where the measure of gauge fields is defined by 
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d A=n 



N N 



=i j= 



(5) 



The action is quadratic with respect to each 
component, which means that we can update each 
component by generating gaussian random num- 
ber in the heat-bath and the Metropolis algo- 
rithm. 

To investigate the scaling behavior of the string 
scale (a'), we study the Wilson loop operator in 
the IIB matrix model. 

The Wilson loop operator and the large iV" 
behavior have been studied with the light-cone 
string field theory of the type IIB superstringj!). 
The Wilson loop operator (w(C)) is defined as, 



w(C) = tr(«(C)), 



(6) 



where C denotes the closed path and v(C) is de- 
fined as v(C) = Up - ■ -Ufj, = Pc exp(i § dak^A^) 
in the bosonic model. The matrices (U^) are con- 
sidered as the unitary matrices, 



Up = cxp(« / dlAp) 



(7) 



In the ordinary lattice gauge field theory, the 
expectation value of the Wilson loop operator 
which spreads a large area behaves as follows, 



w(I, J) ~ exp(— KI x J), 



(8) 



where / and J are the side lengths of the rect- 
angular loop and K denotes the string tension. 
In the same analogy, we study the Wilson loop 
operator of the bosonic model of the IIB matrix 
model. 

From the scaling relation of two-dimensional 
Eguchi-Kawai model, 



g 2 N ~ Constant, 



(9) 



It is expected that the similar scaling relation also 
holds in the IIB matrix model aslHl 



a 1 ~ ,g 2 -/V 7 ~ Constant. 



(10) 



The exponent (7) should be determined dynami- 
cally from the model. For the large N limit, the 
parameter (g 2 N 1 ) must be fixed in the IIB ma- 
trix model in the same manners as the parameter 
{g 2 N) must be fixed in the Eguchi-Kawai model. 

It is considered that the bosonic model is equiv- 
alent to the D > 2 Eguchi-Kawai model in the 
weak coupling limit. Since the U(1) D symmetry 
rotates all the eigenvalues by the same angle, the 
following expansion is valid in the weak coupling 
region, 



e M e H 



(ii) 



where take constant values due to the U(l) D 
symmetry and A^ are small. The bosonic model 
action can be obtained by expanding the action 
of the ten-dimensional Eguchi-Kawai model in 
terms of A^. When the higer order terms of A^ 
can be neglected, we can obtain the area law of 
the Wilson loop operator as 



w(I x J) = <tr(U fl ---U lx ) IXJ > 



(12) 



iA„ 



< tr(e 
exp(-K(I x J)). 



)lx.T > 



In Ref. , the area dependence of the Wilson loop 
operator has been measured and found the area 
law in four-dimension, D = 4. 
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Fig.l The measurement results of the string tension 



We calculate the eigenvalue of Wilson loop op- 
erator numerically in ten-dimension, D = 10, 



w(C) 



1 

=< — trfe 
N V 



i A,, 



)> 



(13) 



We take the loop (C) as the rectangular (I x J) 
where we select any two direction (//, v) in ten- 
dimension. We consider that in the bosonic model 
the isotropy of the ten-dimensional space-time is 
not broken down spontaneously, since the depen- 
dence of the direction of the rectangular is not 
found. 

From the numerical results, the Wilson loop op- 
erator closes to exponential curve with the large 



size area as 



w(I, J) ~ exp(-Jf/ x J), 



(14) 



where K denotes string tension. We plot the re- 
sults in Fig.l. Then we can estimate the string 
scale (a') as 



a' 2 ~ l/K 2 ~ g 2 N imW = Constant. 
3. Summary 



(15) 



obeys the area law. For the scaling behavior of 
the bosonic model of the IIB matrix model, our 
numerical estimation is 



g 2 N xm{i) = ConstanL 



(16) 



Furthermore the numerical result suggests that 
the large N behavior of the square root of the 
string tension approximately equal to the inverse 
of the extent of the space-time HfTof . 



K 



1/2 



N 



-0.54(1) 



R~ 



(17) 



It means that the Planck scale of the theory has 
the same scaling property of the extent of the 
space time in ten-dimension. The numerical re- 
sult is consistent to the suggestion from the string 
field theory on light-cone frame j{| and the 1/D 
expansion ||. 
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Let us summarize the main points made in our 
calculation. We confirm that the Wilson loop 
operator in the ten-dimensional bosonic model 



